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The behavior of homogeneous and disordered systems with a free boundary is 
described on the basis of group theory in the two-loop approximation directly in three- 
dimensional space. The effect of the free boundary on the regime of the bulk critical 
behavior is revealed. It is shown that the boundedness of the system slightly affects 
the regime of the bulk critical behavior in the case of the ordinary transition, whereas 
this effect is more noticeable in the case of the special transition. Surface critical 
phenomena are described for homogeneous and disordered systems, and the critical 
exponents are calculated in the two-loop approximation. It is shown that the effect of 
impurities is insignificant in the special phase transition, whereas it is more noticeable 
in the ordinary phase transition. The derived critical exponents are compared with 
the computer-simulation results. 
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1 INTRODUCTION 

Any systems studied experimentally inevitably have free surfaces, whose effect is usually 
disregarded in description of critical phenomena. However, ordering processes on the free 
surface can proceed at a temperature differing from the temperature characteristic of bulk 
ordering processes, which leads to change in the critical behavior regime. 

There is a temperature range in which surface effects are decisive and are characterized 
by a certain set of the critical exponents. 

The problem of the effect of the boundedness of the system on the critical phenomena 
was first considered phenomenologically by Kaganov and Omelanchuk [T| and in the 
framework of the microscopic approach by Mills [2] and Wolfram et al. [3]. An exact 
solution for the two-dimensional half-plane was obtained by McCoy and Wu [4] by generalizing 
Onsager's work [5] and showed that the phenomenological theory is valid only qualitatively. 
For three-dimensional systems, Binder and Hohenberg [6] showed that the separated 
direction associated with the existence of the free surface makes necessary a separate 
description of critical phenomena on the surface and in the bulk of the system. Binder and 
Hohenberg [6] analyzed two types of the critical behavior that are possible in semibounded 
systems and are attributed to the fact that the spin ordering on the surface occurs earlier 



than in the bulk (surface transition). In addition, Binder and Hohenberg [6] determined 
the critical exponents with the use of the renormalization group approach. The critical 
exponents for semibounded systems with the scalar order parameter were evaluated in the 
framework of ^-expansion in the oneloop approximation by Lubensky and Rubin [7], who 
obtained the value v = l/2+e/12 for the surface critical exponent of the correlation radius 
and the values rj\\ = 2—e/3 and T]± = l—e/6 for the longitudinal and perpendicular critical 
exponents of the correlation function, respectively. In addition, Lubensky and Rubin |8j 
showed that all surface critical exponents for semibounded systems with an n-dimensional 
order parameter are expressed in terms of bulk critical exponents and surface exponent 
rj = (l/2)e(n + 2)/(n + 8) The multicritical point appearing in the simultaneous bulk and 
phase transitions was analyzed in [9, 10] in the framework of the [9l (TOj £ expansion in 
the two-loop approximation. 

The surface and bulk phase transitions in a semibounded system with the n-dimensional 
order parameter were described in [IT] directly in the three-dimensional space in the two- 
loop approximation. The critical exponents calculated in that work are in better agreement 
with previously obtained Monte Carlo simulation results [US, [20l EJJ. In particular, for 
the crossover exponent, computer simulation provides a value of about 0.5, calculations 
directly at D = 3 yield approximately 0.54, whereas the ^-expansion gives a result of 
about 0.68. 

The renormalization group approach to description of slightly disordered systems 
developed in [IH [EH [l_7] directly for three-dimensional systems made it possible to obtain 
the static critical exponents of unbounded systems in the four-loop approximation. The 
effect of the free boundary on the critical behavior of disordered systems was first examined 
in |18| in the framework of the e-expansion. Similar calculations directly in threedimensional 
space in the two-loop approximation were performed in [T9l l20l |2~T] . However, in the 
latter works, the replica procedure was performed with an error, owing to which the 
corresponding asymptotic series for impurity systems are inconsistent with asymptotic 
series for homogeneous systems. 

In all works mentioned above, it is assumed that the presence of the plane free surface 
slightly affects the bulk critical behavior and values obtained for unbounded systems can 
be used for fixed points of the renormalization group transformation when calculating the 
surface critical exponents. However, this statement should be checked and the effect of 
the free boundary on the bulk critical behavior should be estimated. 

In this work, the effect of the plane free boundary on the bulk critical behavior in 
various phase transitions in both homogeneous and disordered systems is examined in 
the two-loop approximation directly in three-dimensional space. In addition, the surface 
critical exponents for the ordinary and special critical exponents for usual and special 
phase transitions in homogeneous and disordered systems with the free boundary are 
calculated with the inclusion of the corrections obtained. 
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2 HAMILTONIAN OF THE SYSTEM 



Let S = S(x) be fluctuations of the scalar order parameter in the half-space V = R+ = 
{x = (f, z)\f e R " 1 ^ > 0}. The system under consideration is bounded by the z = 
plane, which is denoted as dV in what follows. The Hamiltonian of such a model can be 
written in the form 

H = \ I d D x(r + V 2 )S 2 (x) + \ f d D xAr(x)S 2 (x) + u f d D xS\x) (1) 

+ C } [ d D xS\x), 
2 Jdv 

Here, uq is a positive constant; r ~ |T — T c \/T c , where T c is the temperature of the 
bulk phase transition; c ~ \T — T s \/T s , where T s is the temperature of the surface phase 
transition; and At(x) is the random field of impurities such as random temperature. 
Let us pass to Fourier transforms in the f coordinates, 

H = | o °° dz{± J d^qiro + q 2 )S q S- q + d D - 1 ? Ar ! S 5 S_ J (2) 
+u J d D - l qi d D - l q 2 d D - l q z S ql S q2 S q ^ ql _ q2 ^} 



z=0 



For low impurity concentrations, the random field distribution can be treated as 
Gaussian and be specified by the function 

P[Ar] = Aexp[-^- / Ar 2 q d D -\dz\ (3) 

where A is the normalization constant and Sq is a positive constant proportional to the 
concentration of frozen defects of the structure. 

The application of the replica procedure for averaging over the random fields specified 
by the frozen defects of the structure provides the effective Hamiltonian of the system in 
the form 

roo ,1 e m 

H R = dz{- d D ~\{T, + q 2 ) £ S a q S a _ q 
jo z J a=1 

a, 6=1 

m „ 



a=l ' 



+| / d D ~\S q S. 



z=0 



The properties of the initial system can be obtained in the limit m — > for the number 
of replicas. 
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The presence of the free boundary even in homogeneous systems leads to the new 
properties of the model [TT|. The phase diagram of such systems exhibits the disordered 
phase (SD/BD), surface-ordered bulk-disordered phase (SO/BD), and surface-ordered 
bulk-ordered phase (SO/BO). The lines on the phase diagram separating these phases 
determine three kinds of phase transitions. The transitions from SD/BD to SO/BD, from 
SO/BD to SO/BO, and from SD/BD to SO/BO are called the surface, extraordinary, 
and ordinary phase transitions, respectively. The intersection of these three lines of phase 
transitions forms the multicritical point at which the phase transition is called special 
phase transition. The kind of the phase transition is definitely determined by the relation 
between the constants Co and To. In particular, the condition r << Cq << A, where A is 
the inverse lattice parameter and Cq > 0, is satisfied near the ordinary phase transition; 
therefore, the relations r — > and co/r — > oo are valid at the transition point. The 
conditions near the extraordinary phase transition have the form r << Cq << A and 
Co < 0, and r — > and co/t — > — oo along the phase transition line. Near the surface 
phase transition, c << t << A, r > 0; therefore, c — > and c /r — > oo at the phase 
transition point. For the special phase transition, Co — > and r — > 0; however, Co/r — > 0. 

As shown in [7j, the free propagator in the case under consideration has the form 

G(g; z, z>) = i [e-*°'-'l - ^V^')] , (4) 
zkq l Co + fto 

where k = y/q 2 + r . 



3 BULK CRITICAL PHENOMENA 

The behavior of the system in the critical region is determined by the effective charges at 
the stationary point of the renormalization group transformation, which has the following 
form in the case under consideration: 

S (0) = Z l/2 S ^ Tq = b 2 rZ ^ Uq = tf-D uZ ^ Sq = b *-D SZ ^ 

where b is introduced to reduce the quantities to the dimensionless form. 

The Z factors are determined from the requirement of the regularity of the renormalized 
vertex functions, which is expressed in the normalization conditions 

" () r£l*(*0lk»=o = h z 2 r uMlk ^\ k2=0 = b 4 ~ D u, (5) 



dk 2 

7 2 r (4)1 l4-De yr (2,1)| l2-D/2^ 

6 1 8,bulk : |fe2=0 — O 0, Zjl t ,bulk \k a =0 — O T. 

Here T^\ ulk and Tfl ulk are the four-point bulk vertex functions, T^J lk is the two-point bulk 
vertex function, r{ 2 ^) fc is the two-point bulk vertex function with the insert. 

Let us represent the free propagator in the form of the sum of two terms G = Gb + G s 
. The first term 

G b (q;z,z') = ^-e- k ^- z '\ (6) 

ZKq 
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depends only on the difference between the coordinates of two points of the system and, 
after the Fourier transform in the z coordinate, acquires the form 

Gb = — kr, (7) 

To+P 2 

where p = {q,p z ) is the three-dimensional vector. This expression coincides with the 
ordinary propagator for unbounded systems; for this reason, it will be called the bulk 
propagator in what follows. 
The second term 

depends strongly on the position of both points and is determined by the effective surface 
charge cq; this part of the propagator is called the surface propagator. The limiting values 
of the surface part of the propagator are represented in the form 

G a s rd {q; z, z') = 1 e -ko(z+z') ^ 

for the case of the ordinary phase transition for c /r — > oo and 

G r c (q-,z,z') = ^-e- k ^ +z '\ (10) 

for the case of the special phase transition for co/r — ► 0. 

Since the system is bounded only from one side, the regime of the critical behavior is 
determined by the basic volume of the medium that is far from the free surface [6]. When 
considering the points far from the free surface of the system, one can pass to the limit 
z,z' — > oo. In the limiting case under consideration, Gb >> G s . For this reason, in all 
previous works devoted to the critical behavior of semibounded systems, small corrections 
to the bulk critical behavior that appear due to the presence of the free surface are 
disregarded. Nevertheless, the determination of these corrections is an important problem, 
because any systems studied experimentally are bounded. 

Let us write the Callan-Symanzik equation for the bulk vertex functions: 

® d d m d In d r m \ 

Here, the functions 

determine the behavior of the system in the critical region. 

Figure 1 shows the Feynman diagrams only for the vertex functions of the homogeneous 
system, which are necessary for calculation of the corresponding integrals. In these diagrams, 
the lines corresponding to the bulk and surface parts of the propagator are marked by the 
letters "6"and "s respectively. 
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Phc. 1: Feynman diagrams for the vertex function T^ 4 ) 
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The diagrams Jq°\ j[ ^ and contain only "bulk"lines; for this reason, they coincide 
with the corresponding integrals for unbounded systems: 

J*W, J?> = *£, G( 0, = ^. (13) 

The absolute values are given for all other integrals, because their values for the ordinary 
and extraordinary transitions differ only in sign: 

4 1] = ( 2 J)D-i J Q J J d D ~ 1 qdz 1 dz 2 G b (q;z 1 ,z 2 )G s ( y q;z 1 ,z 2 ) = ^, (14) 

| roo roo r 

J ° = (2?r) D - 1 Jo Jo J d°~ 1( ldz 1 dz2G s (q\ z 1 ,z 2 )G s (q; Z\,z 2 ) = — , 

| roo roo roo r r 

J, = J o I J J d D - 1 qd D - 1 pdz 1 dz 2 dz 3 G b (q;z u z 2 ) 

•G s (-q; Z!, z 3 )G b (p; z 2 , z 3 )G b (q - p; z 2 , z 3 ) = 1.421210, 

■G b (~q; Z!, z 3 )G b (p; z 2 , z 3 )G s (q - p; z 2 , z 3 ) = 0.760073, 

, j rOO rOO roo r r 

A = /o \9n-9 / / / / d D ^ 1 qd D ~ 1 pdz 1 dz 2 dz 3 G b (q;z 1 ,z 2 ) 



J ~ J Z " J J > 



(2tt; 



JO JO 



•G b (-q; z 1 , z 3 )G s (p; z 2 , z 3 )G s (q - p; z 2 , z 3 ) = 0.295748, 

j{ 4) = / / / / / d D ~ 1 qd D ~ 1 pdz 1 dz 2 dz 3 G b {q; z u z 2 )G s {-q; z 1 , z 3 ) 
Jo Jo Jo J J 

■G b (p; z 2 , z 3 )G s (q - p; z 2 , z 3 ) = 0.349384, 



, , j rOO roo roo r r 

J i = (2?r) 2D - 2 Jo Jo Jo J J d D ~ 1 qd D ~ 1 pdzidz 2 dz 3 G b (q; zi,z 2 ) 
•G s (-q; z 1 , z 3 )G s (p; z 2 , z 3 )G s {q - p\ z 2 , z 3 ) = 0.174535, 

•G s (-q; Z!, z 3 )G s (p\ z 2 , z 3 )G b (q - p; z 2 , z 3 ) = 0.177920, 

I roo rOO roo r r 

_ / / / / / jD-liD-l„ 



Jl = (2?r) 2D - 2 Jo Jo Jo J J qd Pdzidz 2 dz 3 G s (q;z 1: z 2 ) 

•G s (-q; Z!, z 3 )G s (p; z 2 , z 3 )G s (q - p; z 2 , z 3 ) = 0.112210, 
<#> = -J^o-^f "jf j \d»-\d^pdz,dz 2 dz 3 G b {q + k;z u z 2 ) 



(2tt 

1 d 



•G b (p; zi, z 2 )G s (p + q; z x , z 2 ) fe2=Q = 0.333481, 



Gi 2) = f f / / d^qd^pdz^GM + *; Zl , z 2 ) 



(2tt 

■G s {p- Zl , z 2 )G s (p + q; Zl , z 2 )\ k2=Q = 0.160209, 



Gj 3) = -^L-^A J J d D - 1 qd D - 1 pdz 1 dz 2 dz 3 G s (q + k;z 1 ,z 2 ) 



■G s (p; Zl , z 2 )G s (p + q; Zl , z 2 ) fc2=Q = 0.822467, 
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jspec 
J 

r~iard 

Gspec 


jord 
jspec 



7(0) , 2 r(l) , t(2) 

j -|- zj -t- j , 



G, 



(0) 


,(0) 



3G 



(i) 



.(1) 



3G 



(2) 


(2) 



G 



(3) 

5 



.(3) 



_ 2 _ 2 + jw + 4J w _ 2 jw _ 2 J{^ + J\ 



r(2) 



r(3) 



(4) 



r(5) 



r(6) 



j{ o) + 2 j« + 2 j; 2) + 4 3) + 4 j; 4) 



2J 1 (5) + 2J 1 (6) + J 1 (7) , 



jord^jordy . 

Gord I/ jord\2 
o /Wo J 



0.676982, 
: 0.065781, 



jspec if jspec\2 

J l /Wo J = 
/(^o PeC ) 2 



Gspec 




0.747187, 
= 0.089463. 



Without the surface effects, the ratios of the integrals are 

J/(J ) 2 = 2/3, G/(J ) 2 = 2/27. 



(15) 



Redefining the effective interaction vertices v i = u • Jo and v 2 = 5 ■ J , one arrives at the 
following expressions for the 0— and 7 functions in the two-loop approximation: 



01 



02 



It 



-(4 - D)vi [l - 3617! + 24tj 2 + 1728(2 Ji - 1 - -G)v 2 



9 

2304(2^ - 1 - -G)vi« 2 + 672(2J^ - 1 - -G)^ 
6 3 

-(4 - D)v 2 \l - 24vt + 8v 2 + 576(2 J - 1 - \g)v\ 



1 



1 



- 1152(2 Ji - 1 - -G)v x v 2 + 352(2Jx - 1 G)v 



22 



(4-/J 
192 (2 J - 1 



12u a + 4w 2 + 288 (2 Ji - 1 - ic)^ 2 



;G)vit» 2 + 32(2 Ji 



3 

1 " \g)v, 



Ji 



(4 - D)Q4G{3vl - 3v!V 2 + v 2 2 ) 
G 



G 



J 2 ' 



Ji 
T 2 

The expressions obtained for the /3-functions are asymptotic series and summation 
methods must be used to extract necessary physical information from these series. In this 
work, the following Borel-Leroy transformation generalized to the two-parametric case, 
which provides adequate results for series appearing in the theory of critical phenomena 
1 22 1, is used: 



f(vi,v 2 ) = £ Ci ltia v?v2 = J e'H b F(v 1 t,v 2 t)dt, 
F^v 2 )=gj-f^v{^. 



(16) 

{ix + i 2 + b)\ 

For the analytic continuation of the Borel transform of a function, the following series in 
the auxiliary variable 9 is introduced: 



00 r ■ 

F(v u v 2 ,6) = J20 k j:—- n ^ 



k=0 



Ui U 2 Uj 1 _|_ 42i fc, 



(17) 
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to which the [L/M] Pade approximation is applied at the point 6 = 1. The [2/1] approximants 
with variation of the parameter b are used to calculate the /3-functions in the two-loop 
approximation. As shown in |22| . such variation of b makes it possible to determine the 
range of variation of the vertex functions and to estimate the accuracy of the critical 
exponents obtained. 

The critical behavior regime is completely determined by the stable stationary points 
of the renormalization group transformation, which can be found from the condition that 
the /3-functions are equal to zero: 

AK,^)=0 (i = l,2). (18) 

The requirement of the stability of a fixed point is reduced to the condition that the 
eigenvalues bi of the matrix 

Bhi ^_m^i w=li2) . (19) 

are positive. It is worth noting that the Pade-Leroy summation procedure is possible 
not for any b values and this significantly limits the possibility of applying the method. 
This limitation is associated with the appearance of the poles of the approximants near 
the solutions of the system of Eqs. (TT8l) : for this reason, it is impossible to determine 
the position of the fixed points. In this work, the parameter b varies from to a value 
beginning with which the determination of the stable fixed point becomes impossible. In 
this range, 20 values of the parameter b are taken for which the fixed points are searched. 
Average values with a certain accuracy determined by the spread in the values for various 
b values are taken as the effective charges at the fixed point. 

The stable stationary point of the renormalization group transformation is determined 
by the values (v° rd * = 0.048 ± 0.002 and v%~ d * = for the ordinary phase transition in 
homogeneous systems and by the values vl pec * = 0.066 ± 0.007 and t>2 Pec * = special 
phase transition. These quantities for impurity systems are v ° rd * = 0.067 ± 0.002 and 
v ord* = Q.033 ± 0.003 for the ordinary phase transition and v[ pec * = 0.071 ± 0.001 and 
y spec* _ q_q]_5 _|_ Q.002 for the special phase transition. For comparison, the stationary 
point of the renormalization group transformation in the two-loop approximation is given 
by the values (v* = 0.046 ± 0.002 and v\ = 0) for an unbounded homogeneous medium 
[27] and by the values v{ = 0.067 ± 0.002 and v* = 0.035 ± 0.003) for a medium with 
frozen structure defects [28] . 

The exponent v characterizing an increase in the correlation radius near the critical 
point (R c ~ |T — T c \~ u ) is determined from the relations: 

- ^l+^y 1 = l\l + 6v 1 -v 2 -lU(2J 1 -l-^G-^)v 2 1 



2 V 2 
+ 96(2j; -1- X -G- ~) Vl v 2 - 16(2 Ji - 1 - )p - -)v. 
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The substitution of the effective charges of the fixed points of the special transition yields 
the values u spec = 0.654 ± 0.006 and z/^ c = 0.669 ± 0.003 for homogeneous systems and 
impurity systems, respectively. For the ordinary phase transition, v ord = 0.635 ±0.007 and 
v imp = 0.684 ± 0.002. For comparison, the critical exponents of the system disregarding 
the free boundary for homogeneous systems and impurity systems are v = 0.632 ± 0.004 
and Vimp = 0.685 ± 0.003, respectively. 

The Fisher exponent 77 describing the behavior of the correlation function near a 
critical point in the space of wave vectors (G ~ k 2+v ) is determined in terms of the 
scaling function 7^ as 77 = 7^(1"*, fjj). In the ordinary phase transition in homogeneous 
systems and impurity systems, rf rd = 0.029 ±0.003 and 77°^, = 0.033 ±0.002, respectively; 
in the special phase transition in homogeneous systems and impurity systems, rf vec = 
0.068 ± 0.008 and rj^ p c = 0.069 ± 0.004, respectively. For comparison, the corresponding 
Fisher exponents disregarding surface effects are 77 = 0.030±0.002 and 77^ = 0.035±0.001 
for homogeneous systems and impurity systems, respectively. The other critical exponents 
can be determined from scaling relations. 



4 SPECIAL PHASE TRANSITION 



In addition to the ordinary critical exponents, which are called bulk exponents in what 
follows, the system is characterized by the surface critical exponents on the boundary 
plane. 

Let us introduce the (N + M )-point correlation function 

N M 

G^ M \x,r) = (HS(x t )l[S(r J ,0)) (20) 



i=l 



3=1 



with N points under the surface and M points on the medium surface, and <^...] means 
thermodynamic averaging with the Boltzmann factor exp(—H[S]). The transition to the 
Fourier transforms in the coordinates parallel to the free boundary yields 

N+M 



(21) 



k=i 



M). 



x, r) exp 



N M N M 

- i j2 pAw - i yi PN+j^j] [ n d %i\\\ [ n dr 3 



i=l j=l i=l j=l 

The renormalization group transformations for quantities on the surface have the form 

S^ = ZI /2 S^, co = 6 2 cZn, (22) 



where = S(z) ( 



z=0 



To renormalize the correlation function, the following relation can be written: 

GiT\r,u,6,c) = Z-( n + m ^Z- m/2 G^ m \t ,u ,S ,c ). (23) 
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We introduce the correlation function with inserts in the form 

qi d h 



G (w^) (ro ^ ^ 0)Co) 



G^ M \r ,u ,5 ,c ). 



(24) 



The renormalization group transformations for the correlation function with the inserts 
have the form 

GiT' J ' h) (r, u, 5, c) = Z -^ M ^Z^ M/2 Z^G^ M ^\r Q) u , 5 , c ). (25) 

The expansion of G^ n ' m \t , uq, 5o, cq) in the lowest expansion order in the powers of uq 
has the form 



G^ m \p;t ,u ,5o,c ) 



1 



+ O(u ,S ). 



(26) 



Co + VP 2 + r o 

Similar to bulk vertex functions, the normalization conditions are taken such that the 
surface two-point correlation function for zero external momentum p = coincides with 
its zeroth approximation. Correspondingly, the substitution To — > r, Co — > c provides 

1 



G^ 2 \p;r,u,5,c) 



(27) 



^Gg?(p;T,M,c) 



1 



" " 2^(^+0)2 

For the correlation function with the insert, the normalization condition has the form 

1 



G^\ P ;r,u,5,c) 



(28) 



p=° (c+v^) 2 ' 
The surface Z factors Z\ and Z\\ are determined from the relations 

Z l Z = -2v^(v^ + c) 2 ^G^ (p; r (r, u, 5), « (t, «, 5), 5 (r, u, 5), c (c, r, «)) |^29) 



Z,| = -ZxZ(v^ + cf-^-G™ (0; r (r, it, 5), «o(r, «, 5), 5 (r, «, 5), c ) 

CCq 



CO=Co(c,T,u) 



In order to simplify the calculations, it is convenient to represent these relations in the 
form 



(Z.Z)- 1 = lim 



p^o p dp 

d r 



G (0 ' 2) (p; r (r, u), u (t, u), c (c, r, «)) 



(30) 



Z ll 1 = Z ^ Z TT~ G {0 ' 2 \0;to(t,u,5),Uo(t,u,S),5o(t,u,5),Co) 
1 c*c L 



co=co(c,t,m,<5) 



Let us use the decomposition of the correlation function into the free part 



G (0 < 2) (p;t ,0,0,co) 
and fluctuation corrections cr (p; r , u , S Q , c ), i.e 



= c + k 



G^ 2 \p;t ,u ,S ,c )] 1 =c + k - a (p;T ,u ,5 ,c ). 



(31) 



(32) 
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Phc. 2: Feynman diagrams for the correlation function G^ ' 2 ^ of the homogeneous system. 



To calculate a (p; u Q , So), the total propagator between two points on the surface is 
expressed in terms of the free energy E(p; r , Uq) determining the fluctuation corrections 
to the two-point vertex function as 

G {0 ' 2 \p;t 0j uo,6o,c ) = s \G\ s + s \GTG\ s , (33) 
T = S(l — GE)- 1 . 

The marks s \ and | a show that the left and right points of the propagator are localized on 
the surface. The substitution of Eq. (l32l into Eqs. (l33l ) yields 

(?o(p; r , w , Sq, c ) = 9 °^ 9 ° = g Zg + g EGEg - s \G\ s (g Eg ) 2 + 0(S 3 ), (34) 

1 +s \(j\ s go-L go 

where 



^o(p; ^) = (c + k )G(p; z, 0) = e~ k ° z , k = yV + p 2 . (35) 

The Feynman diagrams for G^ 0,2 ) without impurity vertices in the two-loop approximation 
are presented in Fig. 2 and their analytical expression has the form 

G^ 2 \p; r , 0, 0, c ) = c + fc - E C<(p) - -^T"- (36) 

L J i= i Co + fco 

In this expression, Ci(p) is the contribution from the diagram Bi presented in Fig. 2. In 
this figure, the line marked by D means the propagator 

G D {p, z, z') = 1 \e- k ^ z ' z '\ - e- k ^ z+z 'A , k = Jl+p 2 . (37) 
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The surface vertex function is represented in the form Cq = c + 5c, where .c is the 
shift containing ultraviolet divergences appearing in analysis of surface quantities similar 
to the free energy for bulk quantities. The substitution of Eq. (132)) into (1331 yields 



Z\Z k + c + 5c — a{p) 



p=0 



Vr + c. 



From this relation, 



5c = {[ZiZ]- 1 - l}(v^ + c) + a(0). 
Correspondingly, the correlation function has the form 

G (0 ' 2) = \ko - (y/T + [Z x Z\-\^ + c) - {trip, c + 5c)- (t(0, c + <fc))l ~\ 



(38) 
(39) 

(40) 



This expression provides a scheme for calculating the value of each diagram at p = 0. 
This statement is also valid for subdiagrams. 

For the special phase transition, setting c = in Eq. (39) and expanding 5c in the 
powers of the effective charges, one can obtain the first-order correction 

5c (1) = {[ZxZ]- 1 }^ + a(0). (41) 

The correction obtained is used to determine the first expansion terms for a{p, c + 5c): 



a(p,c + 5 c) 



= a {p, + — 

c=0 <jc 



c=0 



So®. 



(42) 



The Callan-Symanzik differential equation of the renormalization group in the critical 
region has the form 



db +Pl 3 Vl l2 dv 2 2 V+ 2 Vl { ^ Vc > c dc^ ren ' sv ' 



where 



T 



'u— + Ps—)^Z 1 , m = (/3 v —+/3 s —)hxZ ] \. 



du 



35 



du 



05 



(43) 
(44) 



The expression for ri\\ p = rj( spec ^ + r][ spec ^ for disordered system has the form 



■^B 3 - -B + B B - -F 2 ' 



s/ " ' -12B oVl + 12B v 2 + 288(5! + B 2 + ^B 3 - ~B + B B - ~B*)vf 
\ 2\S(Bj B, + -B-, - ^B + B B - \bI)v 2 2 



3 

-576(5! + B 2 + jp 3 + B B - \b 2 ) Vi v 2 . 



(45) 



Here, as well as in the description of critical phenomena, v\ and v 2 are the effective charges 
and 



B - 1 9 B 

Jo up p 

F- 1 92 B 



B,; = — TT^Bj 



^2=0' 



(46) 



J <9p 2 <9c 



c =0,p=0 
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1,2,3). 



The integrals are equal to 

B = 0.5, B 1 + B 2 = 0.323414, B 3 = 0.817579, % = -0.316936. (47) 

The summation of the asymptotic expressions by the Borel-Leroy method and the 
substitution of the effective-charge values at the fixed points of the renormalization 
group transformation provide the surface critical exponents m = —0.239 ± 0.008 and 
V\\imp = —0.279 ±0.009 for homogeneous systems and impurity systems, respectively. For 
comparison, the corresponding critical exponents disregarding the effect of the free surface 
on the bulk critical behavior are r]^ pec) = -0.189 ± 0.001 and = -0.234 ± 0.005 for 

homogeneous systems and impurity systems, respectively. 

The other critical exponents can be found from the scaling relations [TT] : 



v± = 0.5(77 + V\\), Pi = 0.5z/p - 2 + 7711), 7ll = v{\ - 7711), (48) 
7l = p {2 - rj ± ), #1 = — — — , 6 n 



D-2 + 7711 D-2 + 7711 

For homogeneous systems, the critical exponents are 

rfl ec = -0.09 ± 0.01, (3[ pec = 0.249 ± 0.009, ° = 0-81 ± 0.03, (49) 
7 f ec = 1-37 ±0.05, 5^ = 6.5 ±0.8, Sff* = 4.3 ± 0.1. 

For impurity systems, the surface critical exponents acquire the values 

VTZ P = -0-105 ± 0.007, figL = 0.241 ± 0.007, YAtL = 0.86 ± 0.02, (50) 



-imp u.iuj i u.uui, Mlimp u.iti J- u.uui, 

limp 

££ = 1.41 ± 0.04, = 6.8 ±0.4, 5 s ^ mp 



5 ORDINARY PHASE TRANSITION 

The ordinary phase transition requires consideration of the properties of the system in 
the limit c = c/^/r — > 00. However, it is very difficult to examine the function G^ N,M > 
in this limit. As shown in [23l EH [25j [26] , renormalization group analysis is significantly 
simplified for the function 

N M 

G^ M \ x ,r) = (n^)nw))). ( 5i ) 

i=i j=i 

where d n is the normal derivative to the boundary plane. 

Let us denote the correlation function as G^' M ^ = lim-^^ & N ' M \ fluctuation-correction 
function a^ D \p\ r, w , 5q) = cr(p; r, u , 5q, 00), and the free propagator 

G D = Jim Go = -J- \e- kolz - z ' 1 - e -*o(*+*0l . (52) 



2k, 
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Similar to the case of the special phase transition, 

G^ (p; r ,u 0} So) = -k + (p; r, u, 6) , 
Here, a^ D \p; r, uq, 5q) = gT(Go)g, where T(G) is defined in Eqs. (|33~1) . and 



-kz' 



d_ 

dz 



G D (p; z,z') 



z=0 



(53) 



(54) 



The surface renormalizing factor Z 1oo (m, 5), which describes the renormalization group 
transformation of the derivative of fluctuations of the order parameter with respect to the 
normal to the free surface, is introduced as 



(d n S) ren — [Zi i00 Z] 1 / 2 d n S. 
In this case, the renormalization of the correlation function has the form 



(55) 



G { ^({pY,{^t,u } 6) = Z 



-{N+M)/2 y-M/2 



z -M/2 rgwA,) m {z , }) _ (0)1 . (56) 



The normalization conditions are represented as 

d_ 

dp 



GS m (0; r, u, 6) = 0, — Gg% n (p; r, u, 6) 



1 



p 2 =o 2y/r' 

The surface renormalization factor Zi j00 (u,8) is determined from the relation 



(57) 



Ft d 

Zi,oo(u, 8)Z(u, 5) = - lim — 

p-+o p dp 



The scaling function is introduced as 

771,00 



<91nZi j00 dlnZi.c 
Pi ~ r P2- 



dv 9 



(58) 



(59) 



The surface critical exponent rj°f d for the correlation function parallel to the free surface 
is determined from the relation 1 25 1 



(60) 



Tjff* = 2 + 7l lt00 (vl V* 2 ) + T](vl V* 2 ), 

where (v*,v$) is the fixed point of the renormalization group transformation. 
In the two-loop approximation, the scaling function has the form: 



2 - 125 vi - !2B Q v 2 - 288(5i + B 2 



2— 3_ 1— 2n 9 
-B 3 + -So - -B )vl 



-288(5! + B 2 + ^ 3 + ^B - \b\)vI 



(61) 



-576(5! + B 2 



^B 3 - -Bl)vxv 2 , 



Br 



0.410227 B~l + B° 2 =0 B\ 



-0.630378. 
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The substitution of the effective charges at the fixed point yields the critical exponents 
rju rd = 1.51±0.04 and f]°^ p = 1.64±0.06 for homogeneous systems and impurity systems, 
respectively. The other surface critical exponents have the values 

v ^ d = 0.77 ±0.02, [3° rd = 0.79 ± 0.02, 7 °[ d = -0.32 ±0.01, (62) 
7 f d = 0.78 ± 0.03, 5^ d = 1.98 ± 0.06, 5^ d = 0.60 ± 0.02, 
vlLp = 0-84 ± 0.03, fl% d p = 0.91 ± 0.04, 7l TL p = -0.44 ± 0.04, (63) 
-CL = 0.80 ± 0.03, 5^1 = 1.88 ±0.05, 5? d = 0.51 ± 0.03. 



6 CONCLUSIONS 

The calculations reported above reveal the effect of the free boundary on the bulk critical 
behavior of homogeneous systems and disordered systems. The critical exponents v ord = 
0.635 ± 0.007 and rf rd = 0.029 ± 0.003 obtained for the ordinary phase transition in 
homogeneous systems coincide within the errors with the respective critical exponents 
v = 0.632 ± 0.004 and r\ = 0.030 ± 0.002 obtained disregarding the boundary effects. 
Therefore, the presence of the free plane boundary can be disregarded when describing 
bulk critical phenomena near the point of the ordinary phase transition. 

In agreement with expectations, in the special phase transition in homogeneous systems, 
the effect of the free boundary leads to a more noticeable difference of the bulk critical 
exponents v spec = 0.654 ± 0.006 and r) spec = 0.068 ± 0.008 from the respective critical 
exponents for the unbounded systems. However, the difference between indices begins 
with the second place after the decimal point, which is within the errors of experimental 
methods and cannot be detected experimentally. 

A similar pattern is observed for disordered systems. The bulk critical exponents of the 
system with the free boundary coincide with the critical exponents of unbounded systems 
in the ordinary phase transition = 0.684 ± 0.002 and rfi* = 0.033 ± 0.002 for the 
ordinary transition cf. v imv = 0.685 ± 0.003 and rj imp = 0.035 ± 0.001 for the transition 
in the unbounded system). For the special phase transition (z/^ c = 0.669 ± 0.003 and 
Vifnp = 0.069 ± 0.004), the difference in the values of the critical exponents is more 
noticeable; however, it is also within the experimental errors similar to homogeneous 
systems. 

As mentioned above, the difference between the bulk critical exponents is small so 
that it can hardly be determined experimentally. For this reason, the effect of the free 
boundary can be disregarded when calculating the bulk critical exponents. However, the 
effect of the position of the fixed points on the surface critical exponents is much stronger 
and it is necessary to take into account the shifts of these points. 

Table 1 presents the surface critical exponents obtained in this work, as well as the 
surface critical exponents obtained by disregarding the shifts of the fixed points in [TJ] 
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TaSjinna 1: Surface critical exponents for the special transition in homogeneous systems 



A 


7i 




0.249 ± 0.09 
0.257 ±0.006 
0.18(2) 
0.2375(15) 
0.22 


1.37 ±0.05 
1.29 ±0.07 
1.41(14) 
1.328(1) 


This work 
Diehl & Shpot[II] 
Landau& Binder |29j 
Ruge & Wagner [30] 
Vandruscolo h ,a,p.|31j 



and the Monte Carlo computer simulation results for the special phase transition. As 
seen in the table, the results obtained in this work by taking into account that the fixed 
points of the renormalization group transformation are shifted due to the effect of the free 
boundary differ from similar results obtained under the assumption that the boundedness 
of the system does not affect the bulk critical phenomena. However, the ranges of the 
critical exponent values strongly overlap if the errors are taken into account. Comparison 
with the Monte Carlo computer simulation shows that the exponent fix obtained in this 
work is in better agreement with the computer experiment. The exponent 71 is in good 
agreement with the computer simulation results both with and without the effect of the 
boundary on the bulk critical phenomena. 

Comparison shows that the surface critical exponents of the special phase transition 
for homogeneous and disordered systems almost coincide. Differences that are observed 
are likely due to a low order of perturbation theory and will decrease with increasing the 
number of terms retaining in the expansions of f3- and 7-functions. Therefore, the effect 
of point frozen defects of the structure on surface critical phenomena in the special phase 
transition is weak and can hardly be detected experimentally. 

Table 2 presents the surface critical exponents obtained in this work, as well as the 
surface critical exponents obtained disregarding the shift of the fixed points in [TT] and the 
Monte Carlo computer simulation results for the ordinary phase transition. According to 
the tables, the results obtained in this work by taking into account that the fixed points of 
the renormalization group transformation are shifted due to the effect of the free boundary 
are in good agreement with all results of the Monte Carlo computer simulation. It is worth 
noting that the critical exponents of the ordinary surface transition disregarding the free- 
surface-induced shift of the fixed points are inconsistent only with the results reported in 
[29]. Similar to the special phase transition, the effect of the free surface on the position 
of the fixed points cannot likely be detected experimentally. However, in contrast to the 
special transition, the effect of impurities in the ordinary phase transition is stronger and 
can likely be revealed in experiment. 
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Ta6iiHiia 2: Surface critical exponents for the ordin ary transition in homogeneous systems 



Pi 


7i 




0.79 ±0.02 


0.78 ±0.03 


This work 


0.85 ±0.05 


0.74 ±0.02 


Diehl & Shpot[UJ 


0.78(2) 


0.78(6) 


Landau& Binder [29] 


0.807(4) 


0.760(4) 


Ruge & Wagner [30] 


0.79(2) 




Kikuchi & Okabe[32] 


0.80 ±0.01 


0.78 ±0.05 


Pleimling & Selke[33] 
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